A coherent double-quantum-dot electron pump 
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We have investigated the transport characteristics of an electron pump consisting of an asymmetric 
double quantum dot at zero bias voltage which is subject to electromagnetic radiation. Depending on the 
energies of the intermediate states in the pumping cycle, electrons can be transferred through the dots 
incoherently via sequential tunneling or coherently via co-tunncling. The dc transport through the system 
can be controlled by the frequency of the applied radiation. We concentrate on resonant tunneling peaks 
in the pumping current that are interesting to observe experimentally. 
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A quantum dot can be thought of as an artificial atom 
with adjustable parameters. It is of more than funda- 
mental interest to study its properties under various cir- 
cumstances, e.g. by transport experiments M. By con- 
sidering a double-quantum-dot system, the analogy with 
real atoms can be stretched to include artificial molecules. 
The analogue of the covalent bond is then formed by an 
electron which coherently tunnels back and forth between 
the two dots. By applying electromagnetic radiation with 
a frequency equal to the energy difference between the 
time-independent eigenstates of the double-dot system, 
an electron can undergo these so-called spatial Rabi oscil- 
lations even when the tunneling matrix element between 
the dots is small . Recently, several time-dependent 
transport measurements on quantum-dot systems have 
been reported pL most of them aimed to characterize 
the discrete states in the dots. It has also been sug- 
gested to make devices from quantum dots. Examples 
of such applications are pumps that transfer electrons 
one by one by using time-dependent voltages to alternat- 
ingly raise and lower tunnel barrier heights || , or systems 
in which coupled quantum dots (or quantum wells) are 
used for quantum-scale information processing ||. Be- 
cause the transport mechanism in the abovementioned 
pumps is determined by sequential tunneling, electrons 
are pumped incoherently. In this paper we describe an 
electron pump that can also transfer electrons coherently 
by means of a co-tunncling process. In this case, the sys- 
tem switches coherently between two states each time 
an electron is pumped through the dots. The number 
of electrons in the double dot does not change, however, 
resulting in less current noise. Alternatively, this pro- 
cess can be seen as coherent transport through a double- 
quantum-dot qubit. 

The system we consider consists of two weakly coupled 
quantum dots A and B connected to two large reservoirs 
L and R by tunnel barriers (see Fig. |l|). The leads are 
assumed to have a continuous energy spectrum and the 
bias voltage is zero. Each individual quantum dot is in 
its ground state and in the remainder of this paper we 
will only concentrate on the case of transitions between 
ground states of the dots. Moreover, we assume that 
the double-dot system is asymmetric in the ground state 




FIG. 1. Asymmetric double-dot system at zero bias volt- 
age. The dots A and B are coupled to the reservoirs L and R 
by tunnel junctions. An electron in the left dot can tunnel to 
the right by absorbing a photon. 



1 0, 0) , where \n,m) denotes a full many-body state with 
n extra electrons in the left dot and m in the right dot, 
respectively. This asymmetry entails that the energy of 
the state |1, — 1) is much higher than the energy of the 
state | — 1,1). In that case an electron can be excited 
from the left dot into the right one, but the probability 
of the reverse process occurring can be neglected. By 
fabricating two dots of different sizes and adjusting the 
energy levels with a gate voltage, such an asymmetry 
can be easily realized. The energy difference between 
the ground state \a) = |0,0) and the excited state \b) = 
| —1,1) is denoted by eq. It can be tuned by an external 
gate voltage, and we assume it to be much larger than 
the transition matrix element between the dots T, i.e. 
£o ^> T. In this case dc transport through the double 
dot is blocked. 

This situation changes if we apply electromagnetic ra- 
diation to the system. Assume that a time-dependent 
oscillating signal is present on the gate electrode, so that 
the time-dependent energy difference between states \a) 
and 1 6) becomes e(t) = e + ecosojt, where e is the 
amplitude and to the frequency of the externally ap- 
plied signal. When the frequency of the applied radia- 
tion matches the time-independent energy difference be- 
tween states | a) and \b), i.e. if ui « eq, it is possible 
for an electron from the left dot to tunnel to the right 
one. In principle, this electron can now leave the sys- 
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tem by tunneling to the right lead, resulting in state 
|c) = | — 1,0). An electron from the left lead can then 
tunnel to the left dot, thus restoring the ground state. 
Effectively, an electron has now been transferred from 
the left electrode to the right one. This transport cy- 
cle, |0,0) -> -> | - 1,0) -► |0,0), is not the 
only one. Another possible sequence, in which the sys- 
tem passes the intermediate state \d) — |0, 1), is given by 
|0,0) ^|-1,1)^|0,1)^|0,0). 

The details of the transport mechanism of a pump- 
ing cycle depend on the energies of states |c) and \d), 
as shown in Fig. 0. If the energies of these intermediate 
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FIG. 2. Energy diagrams of the double-dot system in which 
(a) only incoherent and (b) only coherent tunneling occurs. 
Indicated are the four possible states of the system and their 
respective energies: \a) = |0,0), |6) = | — 1, 1), |c) = | — 1,0), 
and \d) = |0, 1). 

states are somewhere between those of states \a) and \b), 
the double-dot system relaxes to the ground state via the 
sequential (and thus incoherent) tunneling processes de- 
scribed above. As shown in Fig. || (a), the four tunneling 
processes are described by the rates \ T R \ , and 
F R , respectively. Here, units are used such that h = 1. 
In the opposite case when the energies of the intermedi- 
ate states are higher than that of state \b), Fig. || (b), 
transport via sequential tunneling is blocked. However, 
transport is still possible via inelastic co-tunncling of elec- 
trons [Q. When the system is in state \b), two electrons 
can tunnel simultaneously, one going from the left reser- 
voir to the left dot and the other going from the right dot 
to the right electrode. In this process one of the interme- 
diate states is virtually occupied. The necessary energy is 
provided by relaxing the system to the ground state |a), 
thereby releasing an energy E. In the following, we will 
describe these two different mechanisms quantitatively, 
starting with the sequential tunneling regime. 

We use the density-matrix approach developed in 
Ref. gj. Disregarding tunneling to and from the leads, 
the Hamiltonian of the two-level system is given by: 
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— hi^o + ecoswi) 



T 



|(eo + ecos wt) 



(1) 



where T describes the coupling between the dots. We 
are interested in the particular case where the energy 



level spacing is large, i.e. eg ^> T. Rewriting Eq. ([j]) in 
the basis of eigenstates of the time-independent (e = 0) 
Hamiltonian and using the fact that £<e , we obtain: 



y , -\E Tcosojt 



T cos Lut \E 



(2) 



where E = y £q + AT 2 « eo is the renormalized en- 
ergy level spacing and T = Te/E. From Eq. (g) we see 
that there are small time-dependent matrix elements that 
cause mixing between the two states. Using the equation 
of motion for the density matrix, idp/dt = [H,p], includ- 
ing tunneling to and from the leads, and disregarding 
rapidly oscillating terms, we obtain the equations of mo- 
tion for the density-matrix elements which are valid near 
resonance: 



Paa = +i\T{p ab - pba) + p cc + T^' p dd 

Pbb = -i\T{p ab - Pba) - F T Pbb, 

-n(2) r>( 2 ) 
Pcc = 1 R Pbb - 1 L Pcc> 



(3a) 
(3b) 
(3c) 



Pab = i\T{p aa - pbb) + i(E - uj)p ab - \T T p a b, (3d) 



where T T = T£ + T^> and p aa , p bb , p cc , and p dd = 
1 — Paa — Pbb — Pcc denote the probabilities for an electron 
to be in the states \a), \b), |c) and \d), respectively. In 
the nondiagonal elements p a b = Pb a , only the frequency 
parts pab exp(iwt) and pba exp(— iujt) that contribute to 
the resonant current peak have been retained. 

The average current through the system is given by: 



(2) 



1 R Pbb 



r (l)„ 
1 R Pdd- 



(4) 



In Fig. g we have plotted a resonant current peak as a 
function of the applied radiation frequency for different 
values of T/T, where T = = = = . 
The current peak height initially grows with increasing 
r, its width remaining almost constant. The best reso- 
nance condition occurs for T = T. When T increases be- 
yond this value, the height of the peak rapidly decreases 




FIG. 3. The pumping current as a function of the applied 
frequency for different values of T/T, where we have taken 
= r£ } = ll 2) = = r. Calculated for e = O.leo and 
T = 0.1e . 
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while the width increases. Using an approach developed 
in Ref. S, which is valid near resonance, we are able 
to derive an analytical expression for the shape of the 
current peak in this regime. We find that the peak is a 
Lorentzian: 



In 



I/e = I max w 2 /[w 2 +4(iu-Ej 2 
= T 2 T T /[T 2 T + (2 + a)T\ 



T 2 T + (2 + a)T , 



(5a) 
(5b) 

(5c) 



where a 



/r (i) 

LIB 



T 



(2) /r (2) 



B 



In the limit of small T, 



i.e. in the limit of small radiation amplitude e and small 
overlap T, the height of the current peak is proportional 

2 2 

to T ; Imax — T /Tt, whereas its width is constant in 
this limit; w = Ft- This concludes the discussion of the 
incoherent tunneling regime. 

In the coherent mechanism, two electrons tunnel si- 
multaneously; one going from the left lead to the left 
dot, and one from the right dot to the right lead. Be- 
cause the transport occurs via the virtual occupation of 
a state with a large electrostatic energy, these two tun- 
neling events cannot be treated independently . Using 
Fermi's Golden Rule we obtain the zero-temperature co- 
tunnel rate 101: 



ct 



2C [ln(l + E/A L ) + ln(l + E/Aj 
A L + A R + E 
CE CE 



A L (A L +E) A R (A R + E)' 



(G) 



where C — 2ir\tLtR\ 2 i'Li / R, and where the matrix ele- 
ments for tunneling trough the left and right barier, tL,R, 
and the density of states in the left and right electrode, 
Vl,Bi are assumed to be energy independent. In the limit 
of large charging energies, Al, Ar 3> E, which we will 
consider from now on, Eq. (|^) reduces to a simple expres- 
sion: 



2%\t L t R \ v L v R 
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The prefactor C is proportional to the product of the 
individual tunnel rates through the left and right bar- 
rier, and the energy denominators reflect the fact that 
the tunneling occurs via the virtual occupation of two 
states. Because co-tunneling is a second-order process, 
this co-tunnel rate will be much smaller than the sequen- 
tial tunnel rates of the incoherent regime. In the absence 
of radiation, the average current through the system is 
given by: 



I/e = T ct p bb , 



(8) 



taking into account the radiation-induced co-tunneling 
processes. In contrast to the situation in the incoherent 
mechanism, these density matrix elements are the only 
diagonal ones because the states |c) and \d) are now oc- 
cupied only virtually. 

As before, we choose the eigenstates of the time- 
independent Hamiltonian as the basis for our calcula- 
tions. The time-dependent Schrodinger equation then 
becomes: 



idi) a (t)/dt = H(t)tp a (t), 



(9) 



where H. is given by Eq. (g) and a = a,b. Expanding the 
wave functions ipa(t) in harmonics we obtain: 



—iE a t \ ^ ~>n — incut 



(10) 



In the case of a small radiation amplitude, e <C Eq, and 
near resonance, u> ~ E, it is a good approximation to 
retain only the n = and n — 1 terms for ip a (t), and the 
ri = and n = — 1 terms for ipb(t)- This approximation 
entails that we only take into account single-photon tran- 
sitions between the states \a) and \b). The wave functions 
then become: 



JE a t 



where the coefficients £ and £ are given by: 



£ = (6E + W)/y/2W(W + 6E), 
( = T/y / 2W(W + 5E), 



(11a) 
(lib) 



(12a) 
(12b) 



where SE = E - to and W 2 = SE 2 H 
when lo — E, the coefficients are £ 



-T 2 



energy E a in Eqs. (11a) and ( lib ) is given by 
E a - -E b 



At resonance, 
C = hV2. The 
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(13) 



(7) reflecting the fact that the energy levels —\E and ^E- 



slightly repel each other when u> « E. 

With the aid of the closed-time-path Green function 
technique [^o| , which was applied to quantu m-do t sys- 
tems in Ref. fill, the wave functions Eqs. (11a) and 



(lit) are used to obtain the equations of motion for the 
density-matrix elements near resonance. Without going 
into the details of this calculation (they can be found 
in Ref. pi), we simply present here the resulting set of 
equations: 

Pbb = -i\T{p ah - p ba ) - T ct pbb + C 2 C 2 r c t, (14a) 



Pab 



i\T(p aa - p bb ) + i(E - uj)pab - h^ctPab, (14b) 



which is zero, since p bb = 0. In order to calculate the 
current in the presence of radiation, we therefore need 
equations for the density matrix elements p aa and p bb 



where p aa + p bb = 1, p ba = p* ab , and where we have used 
u> w E and E a w —E/2 to calculate the dissipative terms. 
Clearly these equations are very similar to Eqs. (^) and 
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(pq). The crucial difference with the equations for the 
density matrix in the incoherent mechanism, however, is 
that the extra tunnel rate in Eq. (14a) depends on the 
applied frequency and radiation amplitude via £ and 
It is a Lorentzian centered around ui = E having width 
2T. 

We solve for the stationary solution of these equa- 
tions to calculate the average current through the system 
which is given by the probability pbb to be in the excited 
state | b) times the decay rate of that state: 



I/e 



c 4 )r ct pbb- 



(15) 



In Fig. H the scaled pumping current is plotted for differ- 
ent values of T ct /T. With increasing co-tunnel rate, the 




FIG. 4. The scaled pumping current in the coherent regime 
as a function of the applied frequency for different values of 
r c t/T. Calculated for e = O.leo and T = 0.1s . 

scaled peak height first drops rapidly and then decreases 
slowly to its minimum value, while its width changes only 
slightly. The peak shape is a Lorentzian of the form 
Eq. ( |5a| ) with height I ma x and width w: 



In 



Hi 



(16a) 



TV6(2T +T 2 ct )(AT +T 2 ct )/(8T + (16b) 



ir ct (4T +T 2 ct )/(2T +T 2 ct ), 



reduces to |r ct 



for small T r + and to 



The height I„ 

^r ct for large r ct . A comparison with Fig. || clearly 
shows that the co-tunnel current peak is much narrower 
than the incoherent one. Another distinct feature of this 
co-tunnel peak is the fact that its width changes non- 
monotonically with increasing co-tunnel rate r ct . For 
small values of T ct the width is equal to w — \/6 T « 
2.45 T . It then increases rapidly to reach a maximum of 

2 3/4 55 
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\V2 T « 3.52 T for T c 



1.68 T. 



2 v - j- ~ ^-u- ^ ct 

On increasing the co-tunnel rate further, the width de- 
creases again to v6 T. The width is thus comparatively 
insensitive to changes in the co-tunnel rate. Thi s is due 
to the fact that the extra co-tunnel rate in Eq. (14a) it- 
self has an intrinsic width w = 2T, thus restricting the 
width of the current peak dramatically. 



In conclusion, we considered an electron pump con- 
sisting of a double quantum dot subject to radiation. An 
incoherent and a coherent pumping mechanism has been 
discussed. By deriving equations of motion for the den- 
sity matrix elements of the double-dot system, we cal- 
culated the pumping current in both regimes. Whereas 
the tunnel rates as a function of the external frequency 
are constants in the incoherent mechanism, they are 
Lorentzians in the co-tunneling regime. In both cases 
the current peak is a Lorentzian, but in the coherent case 
the peak width is much smaller than in the sequential- 
tunneling regime and changes nonmonotonically as a 
function of the co-tunnel rate. Experimental realization 
of this device would allow for a systematic study of co- 
herent transport through a solid-state qubit. 
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